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C^ ■ Abstract 

A Hubbard model with a A'^z-fold degenerate single-particle ground state has ferromagnetic ground 
'^ . states if the number of electrons is less or equal to Nj. It is shown rigorously that the local stability of 

C ' ferromagnetism in such a model implies global stability: The model has only ferromagnetic ground 

states, if there are no single spin-flip ground states. If the number of electrons is equal to Nd, it is well 
known that the ferromagnetic ground state is unique if and only if the single-particle density matrix 
is irreducible. We present a simplified proof for this result. 
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Qv ' The problem of ferromagnetism in itinerant electron systems has a long history. It is clear that ferromag- 

netism (as any other ordering in itinerant electron systems) occurs due to the interaction of the electrons, 
or, to be more precise, due to a subtle interplay between the kinetic motion of the electrons and the 
interaction. 1963, Hubbard [|I]], Kanamori [^, and Gutzwiller [^] formulated and studied a simple tight- 

'-^ ■ binding model of electrons with an on-site Coulomb repulsion of strength U. This model is usually called 

C , the Hubbard model. Although the assumption, that a realistic system can be described by a purely local 

repulsion of the electrons is artifical, the Hubbard model became a paradigm for the study of correlated 
electron systems. The reason is that already a pure on-site interaction can produce many ordering effects 
that have been observed in electronic systems. The mechanisms that are responsible for some long range 
order in the ground state of the Hubbard model are probably the same in more complicated (and more 
C^ ■ realistic) models. From a theoretical point of view the Hubbard model is very interesting, because it offers 

the possibility to derive ordering phenomena in a simple model that does not contain special interactions 
favouring this order. 

In this paper we present a result on ferromagnetism in the Hubbard model. This is an old problem, 
which has been studied extensively using various approximative methods. The most simple approach is 
the Hartree-Fock approximation. It yields the Stoner criterion Upp > 1 for the occurence of ferromag- 
netism in the Hubbard model, pf is the density of states at the fermi energy. It is well known that this 
criterion overestimates the occurence of ferromagnetism. There are situations where pfU is infinite and 
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the ground state of the Hubbard model is not ferromagnetic. Ferromagnetism is not a universal prop- 
erty of the Hubbard model. As far as we know it occurs on special lattices and in special regions of the 
parameter space. 

In the discussion of ferromagnetism in correlated electron systems, more realistic models with e.g. an 
additional ferromagnetic interaction between the elctrons [Q] or a Hund's coupling between several bands 
in a multi-band system [||] have been discussed as well. It is clear that in a realistic description of itinerant 
ferromagnetism such additional interactions are present and may favour the occurence of ferromagnetism. 
But it is a challenging problem to derive conditions for the occurence of ferromagnetism in a Hubbard 
model, which does not explicitely contain such interactions. The hope is that results for this model yield 
a main contribution to the understanding of ferromagnetism in more realistic models. 

The first rigorous results on ferromagnetism in the Hubbard model is the so called Theorem of 
Nagaokal^. On a large class of lattices, the Hubbard model has a ferromagnetic ground state if the 
Coulomb repulsion is infinite and if there is one electron less than lattice sites. A very general proof of 
this theorem has been given by Tasaki [^. 

A second class of systems, for which the existence of ferromagnetic ground states has been shown 
rigorously, are the so called flat band models. In 1989, Lieb [||] proved an important theorem on the 
Hubbard model: At half filling and on a bibartite lattice (one electron per lattice site) the ground state of 
the Hubbard model is unique up to the usual spin degeneracy. The spin of the ground state is given by 
5=^||A| — |B||, where \A \ and \B \ are the numbers of lattice sites of the two sublattices of the bipartite 
lattice. When this quantity is extensive, the system is ferrimagnetic. In that case, the model has strongly 
degenerate single particle eigenstates at the Fermi level, Pf is infinite. For a Hubbard model on a transla- 
tional invariant lattice such a model has several bands, one of these bands is flat. Later, it has been shown 
that a multiband Hubbard model for which the lowest band is flat shows ferromagnetism [P] - [|T4|] . These 
lattices are not bipartite, generally they contain triangles or next-nearest-neighbor hoppings. A typical 



example is the Hubbard model on the kagome lattice [11]. 



There are several extensions of the flat band ferromagnetism. The most important result has been 



derived by Tasaki [|15||. He discussed the question, whether the flat band ferromagnet is stable with 
respect to small perturbations. He showed under very general asumptions that for a class of multi-band 
Hubbard models with a nearly flat lowest band the ferromagnetic state is stable with respect to single spin 
flips if the Coulomb repulsion U is sufficiently large and if the nearly flat band is half filled. This local 
stability of the ferromagnetic state suggests its global stability. The class of models, for which Tasaki was 
able to proof this important result consists of models, for which the nearly flat, lowest band is separated 
from the rest of the spectrum by a sufficiently large gap. Therefore one would expect that these models 
describe an insulating ferromagnet. This a general problem for the flat band ferromagnetism as well. The 
flat band models show ferromagnetism if the flat band is half filled or less than half filled. Even if the flat 
band is less than half filled, or if the model has no gap between the flat band the other bands (this is the 
case for the kagome lattice), the system may be an insulator. The reason is that for an entire flat band, 
the systme can be described by localized states as well. Furthermore, the existence of a basis of localized 
states was an essential part of the proofs. 



Another extension of the flat band ferromagnetism are models with a partially flat band. In [ |13[ ] a 
general necessary and sufficient condition for the uniqueness of feiTomagnetic ground states has been 
derived for a model with a degenerate single particle ground state. This result holds only if the number of 
electrons is equal to the number of degenerate single particle ground states. But it does not require a gap 
in the spectrum or an entirely flat band. A partially flat band is sufficient. A Hubbard model with a single 
band far away from half filling is expected to be a conductor. This remains true if the band is partially 
flat. Therefore these models may describe a metallic ferromagnet. A generalization of this result to a 
situation where the number of electrons is less than the number of single particle states has recently been 



published [16]. The main result of that letter is that in a single band Hubbard model with a degenerate 



single partical ground state local stability of ferromagnetism implies global stability, if the number of 
electrons is less or equal to the number of degenerate single particle groud states. Stability is meant here 
in the sence of absolute stability: The ferromagnetic ground state is the only ground state of the system. 

That local stability of ferromagnetism implies global stability has often been assumed but is by no 
means guarantied. It would be certainly useful to know, in which situatons this is the case. The aim of 
the present paper is to generalize the result in to a general Hubbard model with degenerate single particle 
ground states. It is not necessary to have a single band model. It is even not necessary to have translational 
invariance, although this would be a natural asumption. 

Let us mention that Hubbard models with a partially flat band are not only an academic toy model. 
Very recently Arita et al [ p7| ] used such a model to explain the negative magnetoresistance of certain 
organic conductors. They mention that standard band-structure calculations for these materials yield a 
partially flat band. 

The paper is organized as follows. The next section contains the main definitions and results. The 



proof of the result combines the ingredients of the proofs in [ ]13[ ] and. The main part of the proof is the 
choice of a suitable basis. This choice is discussed in Section 3. Section 4 contains a new proof of the 



result in [13]. The proof in [13] used an induction in the number of degenerate single particle states and 



was not intuitive at all. On the other hand, the basic idea why the condition in [13] should be true is 
simple. The new proof is based on this basic idea and is much easier. Furthermore it can be generalized 
to situations where the number of electrons is less than the number of degenerate single particle states. 
This generalization is presented in Section 5. 



2 Main result 



We consider a Hubbard model on a finite lattice with N, sites. The Hamiltonian is 
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X and y are lattice sites. As usual cJo and c^a are the creation and the annihilation operators of an electron 



on the site x with a spin a =|, J,. They satisfy the anticommutation relations 

{Cx<3iCyz\^ 



C. va , Cyx 



5vv5cTt, and 






= 0. The number operator is defined as Uxa = CxaCxa- The hopping matrix T = 

{txy) is real symmetric, and the on-site Coulomb repulsion Ux is positive. We do not need to assume any 
kind of translational symmetry, therefore the lattice is simply a collection of sites. We allow that the local 
Coulomb repulsion depends on x. The total number of electrons is Ne = LAeA(".vT ~^^xi)- 
The Hubbard model has a SU{2) spin symmetry, it commutes with the spin operators 



where ']f is the vector of Pauh matrices, 
/O 1 \ 
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We denote by S{S + 1 ) the eigenvalue of S . 

The eigenstates of the hopping matrix T are cp', the corresponding eigenvalues are £,, / = 1, . . . ,Ns. 
Let £,- < Ej for / < j. In the following we will discuss a Hubbard model with N^ degenerate single particle 
ground states. The energy scale is chosen such that £,■ = for / < Nd- 

The main result of the present paper can now be formulated. 

Theorem - In a Hubbard model with Nd degenerate single particle ground states and Ng < Nd 
electrons, local stability of ferromagnetism implies global stability: The model has only 
ferromagnetic ground states with a spin S = ^, if and only if there are no single spin-flip 
ground states (ground states with a spin 5=^ — 1). 

This theorem is true for any positive Coulomb repulsion Ux- The existence of feromagnetic ground 
states is indeed trivial. Any multi particle state that contains only electron with spin up in single particle 
states cp', i < Nd is a ground state of the Hamiltonian. It is even a ground state of the kinetic part and 
of the interaction part of the Hamiltonian separately. The problem is to show that there are no further 
ground states. The above theorem yields a necessary and sufficient condition for the existence of non- 
ferromagnetic ground states. 



As already mentioned in [|13|], one can use degenerate perturbation theory to generalize the result to a 
situation where the flat part of the band lies not at the bottom of the single-particle spectrum. A concrete 
model, where such a situation occurs, has been investigated by Arita et al [|l^]. They investigated a special 
type of one-dimensional Hubbard model used for the description of atomic quantum wires. These models 
have a flat band that is not situated at the bottom of the spectrum. 

3 Choice of the basis 

The main part of the proof of the theorem is the choice of an appropriate basis. It turns out that the choice 
of the single particle basis used in [ [l3| ] is useful. In this section we give a detailed construction of such a 
basis. Starting point is the representation 

ToC To ) ^^^ 

of the hopping matrix. Here Tq is a positive {N^ —Nd) x (Ns — A'^)-matrix, rankJo = N^ — Nd- C is a 
{Ns —Nd) X Ni-vaairix. This representation of T can be obtained as follows: Since rankr = N^ —Nd, one 
can find N^ — Nd rows (or colums) of T wich are linear independent. We label the con^esponding sites by 
X = Nd + l,... ,Ns. To is the submatrix {txy)x,ye{Nd+i,...,N,}- Since T is non-negative, To is positive. The 
matrix C is given by (ro)^'roi where Tqi = {txy)xe{Nj+i,...,Ns},ye{i,...,Na}- The other matrix elements of T 
are fixed, since T is symmetric and since the other rows of T are linear dependent. 

By construction, the single particle ground states obey r\|/ = 0. This holds if and only if 

A basis of single particle ground states can be obtained by choosing an arbitrary set of Nd linear indepen- 
dent vectors \j/. We choose the basis ® = {\|/; : \j/;(x) = 5xj}. Since txy are real, \|/,(x) are real. This basis 
is not orthonormal. The matrix B = {bij)ij=i^,,,^Nj with bij = LxV,(x)\|/y(x) is positive and the dual basis 
is formed by \|/f (jc) = ^y(B^^),y\|/y(jc). One has Lv¥f (■^)¥i(-^) = ^tj- We introduce creation operators 
for electrons in the state ^i{x) 

«L = l¥.-W4 (8) 



and the corresponding dual operators 

aia = Y.V!{x)cl (9) 

They obey the commutation relations [a,o,aU = 5,j5oj- These creation and annihilation operators can 
now be used to construct multi-particle states. The unique ferromagnetic ground state with Ne = Nd 
electrons and S = St,= Nci/2 is 

¥of = n4lO) (10) 

i 

A general ground state of the kinetic part of the Hamiltonian is given by 

i^"-'"(a)=5l''"(a)^^0F (11) 

where 

S":'"(a)= I ay....i„,;/,...4,n«Mn«4T (12) 

j[...j„,;ii...i„ k k 

This state has Ne = N^ — n + m electrons. In the following I assume that Ng < Nd, i.e. m < n. \|/"''"(a) is 
a state with 53 = (A'^^; -n-m)/2 = Ne/2 - m. It obeys 5+\|/"'™ (a) = if and only if 

Y.^Kh-h,-i^h-in-i =0- (13) 

k 

In that case it is a state with a spin 5 = ^3. We want to derive a condition for \|/"''"(a) to be a ground state 
of the Hamitltonian. A necessaiy and sufficient condition is 

C4C4^'"'(a)=0 (14) 

for all X. For x < Nd one obtains a;|a,|\|/"''" (a) = 0. Therefore OCyi ..y„,;,-, ...,„ = if {71 . . . j,„} is not a subset 
of {/i . . . /„}. It turns out that this fact is important since it simplifies the proof substantially. 

4 The case Ne = N^ 



This case has already been discussed in [|13|]. In the following we obtain a simplified proof of the result 
in [13]. Let us first discuss the stability with respect to single spin-flips. A ferromagnetic ground state is 
called stable with respect to a single spin flip, if there is no single spin flip state with the same energy. I 
derive a necessary and sufficient condition for \|/of to be stable with respect to a single spin flip. A general 
state with a single spin flip can be written in the form 

¥ = £a;;Ky|C^TVOF (15) 

If and only if aj-^t <=< 5yjt, \|/ is the unique ferromagnetic ground state with S = Nd/2, S^ = Nd/2 — 1. If and 
only if T^jOijj = 0, \|/ is a state with S = S^ = Nd/2 — I. Therefore I assume Y^j^j-.j =0. \|/ is a ground 
state if and only if Cx^Cxi^ = 0. 

Cx]Cxl^ = Y^ aj-k^j{x)M^i (x)q|Q|\|/0F (16) 

j,u 
The right hand side vanishes if and only if 

Y,Mfj{x){aj-u^i{x) - aj-iMfk{x)) =0 \/k,l (17) 

i 



I introduce 

%W=La7;'t¥;W (18) 

j 

The condition for aj,k yields 

^k{x)^\fl{x) - ^l{x)^\fk{x) = V^, I, x (19) 

A trivial solution is \jifyt(x) = \|//t(x). It corresponds to ayjt oc 5j^k and has been excluded above. Multiplying 
the condition for ^(x) with Yi{y)Vkiz) ^'^'1 summing over k and / yields 

Py,xPx,z - Py,xPx,z = (20) 

where Py^ = T,j^j{x)\^j{y), Pj,,x = llj'^j{x)wj{y)- If the matrix p^^ ^ is irreducible, the only solution is 
P)',.v = P)',.v- It corresponds to ttyi oc bj^k and has been excluded above. If the matrix P;i ^ is reducible, 
the equation for p^jc has another non-trivial solution. From the non-trivial solution for Pj, v one obtains a 
solution for ttyi from which one can construct easily a solution with Y,j Gijj = 0. Thus \|/of is stable with 
respect to a single spin flip if and only if p^y is irreducible. This is the condition derived previously in 
[p^]. To derive this condition, it was not necessary to use the special single particle basis introduced in 
Section 3. The use of this basis is useful for the investigation of multi spin-flip states. 
Let us now consider a multi spin-flip state \|/"''"(a). It is a ground state if and only if 

£(-!)''£¥;, W¥^P(„^„a;,...7„;i.,i)...%„) = (21) 

Since ay,...y„;,,...,„ is antisymmetric in the last n indices, it is sufficient to sum over all cyclic permutations. 

n+l 

L (-l)'"'L¥7, ix)^\fkXx)aj,...j„-Mi-kn^iM-kr-i = (22) 

Since (Xj^,,,j^-i^,,j^ ^ only if {71 . ..j„} = {/i . ../„}, we obtain for « = 1 (the single spin-flip case) 

\|/^(x)\|/,t' {x){at,k - CLk'-k') = (23) 



With Vfk{x) = Uk,k^k{x) this yields the original condition (19). This means, that with this choice of the 
basis, the functions \|//t(x) are either equal to "^k{x) or vanish. A solution exists, if the set {\|/i(x),^ = 
\ ...Nd} decays in two subsets such that \|/yt(x)\|/^/(x) = if the two factors are out of different subsets. 
This is equivalent to the above condition on the single particle density matrix Px^y For « > 1 we now use 
the fact that the set { ji , . . . , 7,,} is a subset of {^1 , ... , fc„+i }. I choose 72 = ^1 , 73 = ^2, etc in (p2|). With 
this choice only the terms with r >n in the sum over r do not vanish. For r = n the only non-vanishing 
contribution in the sum over ji is ji = /c„+i. For r = n+l one has j'l = /:„. This finally yields 

¥^„(-^)¥i„+, {x){ak^...k„-M-k„ - a^„+ii,...i„-,;i„+,M...^„-i) = (24) 

For some fixed ^i,. . . ^„_i I let at^k = otytyti.../t„_i;/t/ti...yt„_i- The indices ki ... k„^\ are chosen such that (Xk-^k 
does not vanish identically, which is possible since (Xy(...y,,;,,..j„ does not vanish identically. This shows 
that the existence of a multi spin-flip ground state implies the existence of a single spin-flip ground states. 
Therefore the ferromagnetic ground state of the Hubbard model with N^ = Nd electrons is the unique 
ground state (up to the spin degeneracy due to the SU{2) symmetry) if and only if p^y is irreducible. 



5 The case Ne < Nd 

It is now very easy to generalize the second part of the above derivation to the case Ng < Nd- We will 
show that the existence of a multi spin-filp ground state implies the existence of a single spin-flip ground 
state. Let \|/^'"(a), « > 1 be a single spin-flip ground state for Ne=Nd — n+\ electrons. The condition, 
that this is a ground state, yields 

n+l 

£(_1)«'- £ Vy,W¥-t,.Wa7,A.+,.../t„+,-t,..i>_i=0 (25) 

r=l j,e{ku-k„+i}\{k,} 

The sum over j'l is restricted to the set {k\,. . .A:„+i} \ {kr} since otherwise o~.j^-^kr+i...k„+iki...kr-i vanishes. 
The similar condition for a multi spin-flip ground state \|/'"'"+'"^^(a) is 

n+m 

£ (_l)(«+-i)^ £ ^.^ (x)v,,.(x)ay,..,„,;,,.,,..,„,,„,,..,,._. = (26) 

'•=1 jl&{kl,-k„+m}\{k,-,i2-im} 

I let jr = kn+r, r>2. Then the sum over r runs from 1 to « + 1 and the sum over j'l runs over all elements 
of {k{ . . .kn^y] \ {kr}. all other terms vanish identically. One obtains 

n+\ 

£(_!)'"• £ V;iW¥*:,.WcX;ii„^2...«:„^,„;«:,.^l...^„^li:i...i:,._li„^2...«:„^,„ =0 (27) 

'■=1 ;ie{«:i,...«:„+i}\{<:J 

Therefore we can choose a;,;yt,^,...i„^,yt,...i,_, = a;,fc„+2...yt„+„;/t,+,...4„+i,t,...i,_,fc„+2...i„+,„ for some fixed /c„+2, ■■■kn^ 
such that ay[;i,.^,..jt„^,ii..jt,._idoes not vanish identically. This is possible since (X^i...;,„;/ti.../t„ does not vanish 
identically. The corresponding single spin-flip state \|/^'"(a) is thus a ground state iov Ng = Nd — n + \ 
electrons. 



The proof presented here is considerably simpler than the proof in [|13j]. Compared to the proof 
in [ [r6| | it has the advantage that the existence of the single-spin flip ground state is trivial, whereas in 
[ p^ ] a lengthy calculation (hidden in fotenote 13) was necessary to show that. On the other hand, for a 
translationally invariant multi -band system the basis used here is clearly artifical. But if one uses a natural 
basis of Bloch states, it is very hard to construct the single spin-flip states from multi spin-flip states. 

6 Summary and Outlook 

In this paper it is shown that for a general Hubbard model with degenerate single particle ground states 
local stability of ferromagnetism implies global stability of ferromagnetism. To be more precise: If there 
are no single spin-flip ground states, all ground states have the maximal spin. This result holds if the 
number of electrons is less than or equal to the number of degenerate single particle ground states. A 
similar result has been proven for a single band Hubbard model in [|l^] and the present result can be seen 
as a generalization. It holds for a Hubbard model with more than one band, it holds even if the model 
does not have translational invariance. Furthermore, the proof presented here is much simpler than the 



proof in Qlq ] 



The result is important since in many cases it is much simpler to show local stability of ferromag- 



netism than global stability. In [15] showed that under very general conditions the flat band ferromag- 
netism can be extended to situations where the lowest band is not flat but has a weak dispersion. He was 
able to prove that in such a situation the ferromagnetic ground state is locally stable. It would be very 
interesting to obtain conditions under which in that case the ferromagnetic ground state is globally stable, 
i.e. where it is the real ground state of the system. This is clearly a very difficult project. The present 
theorem does not apply since Tasakis model does not have degenerate single particle ground states. But 



one may hope that a generalization is possible. If in a situation with degenerate single particle ground 
states the ferromagnetic ground state is the only one, it is possible that ferromagnetism is stable with 
respect to small perturbations of the Hamiltonian. The main problem is clearly that for a general model 



there is no gap in the single particle spectrum like in the models discussed by Tasaki []15|]. 

From a physics point of view ferromagnetism for models with a partially flat band, as studied in this 
paper, differs from the flat-band ferromagnetism, since a flat-band ferromagnet is typically an insulator, 
whereas models with a partially flat band describe metals. Therefore our new approach is a step towards 
the understanding of metallic ferromagnetism in the Hubbard model. 
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